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Q^ ■ We consider models of light super-weakly interacting cold dark matter, with C(10 — 100) 

C^ ■ keV mass, focusing on bosonic candidates such as pseudoscalars and vectors. We analyze 
the cosmological abundance, the 7-background created by particle decays, the impact on 

1^^ , stellar processes due to cooling, and the direct detection capabilities in order to identify 

O \ classes of models that pass all the constraints. In certain models, variants of photoelectric 

59 ■ (or axioelectric) absorption of dark matter in direct-detection experiments can provide a 

:• ■ sensitivity to the superweak couplings to the Standard Model which is superior to all existing 

• 1— I ! indirect constraints. In all models studied, the annual modulation of the direct-detection 

r> \ signal is at the currently unobservable level of O(10~^). 
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1 Introduction 

The evidence for the existence of non-baryonic dark matter (DM) now comes from many 
sources and ranges over many distance scales pj, from the rotation curves of galaxies, the 
dynamics of clusters, lensing data and the characteristics of large-scale structure, to the 
features of the cosmic microwave background (CMB) fluctuation spectrum and the success 
of big bang nucleosynthesis (BBN). All of these pieces of astronomical data point to a similar 
cosmological density of dark matter, several times that of visible baryonic matter. However, 
this data only probes the gravitational interaction of dark matter and, while it presents 
us with one of the most compelling arguments for physics beyond the Standard Model, 
gaining insight into its non-gravitational interactions remains a primary experimental focus, 
both through underground detectors, particle colliders, and the observation of photon and 
neutrino fluxes from overdense regions in the galaxy and beyond. 

This lack of any direct information on how dark matter may couple to the Standard Model 
(SM) means we are forced to rely on various theoretical expectations. In particular, the 
many successes of standard cosmology motivate a simple thermal mechanism for populating 
the universe with dark matter with a well-defined freeze-out abundance as the universe 
expands. This in turn requires a specific annihilation cross-section which is necessarily non- 
gravitational in origin. The fact that a weakly interacting particle with a weak-scale mass 
has an annihilation rate in the right ballpark pi, T], combined with our expectations for new 
physics at the electroweak scale, has rightly led to the prevailing WIMP paradigm for cold 
dark matter. Nonetheless, persistent problems in understanding the small-scale gravitational 
clustering properties in cold dark matter simulations and on galactic scales has motivated 
variants of this picture where the dark matter may be somewhat lighter, with masses down 
to the keV range. Masses in this range imply a super-weak interaction strength between dark 
matter and the SM sector, indeed many orders of mangitude below weak-scale cross sections. 
This follows from the necessity to have early thermal decoupling of the DM sector, prior to 
the electroweak epoch at T ~ 100 GeV, in order to satisfy the conflicting requirements of 
not having too much energy density in dark matter, and the strong lower bounds on itidm 
coming from the analysis of structure formation. If early decoupling can be achieved, then 
masses in the keV range may withstand these combined constraints and at the same time 
provide a rather attractive mechanism for ensuring the correct dark matter energy density. 

An important feature of keV-scale dark matter is that, unlike the majority of electroweak 
scale WIMP models, it need not be stable against decays to light SM degrees of freedom, 
e.g. photons and neutrinos. Given the super- weak strength of its interaction such decays 
may be strongly inhibited, but nevertheless keV dark matter can be emitted or absorbed 
in astrophysical environments and in terrestrial experiments, and thus is possibly subject 
to additional constraints. The main question to be studied in this paper is whether the 
application of cosmological and astrophysical constraints leaves any realistic chance for the 
direct detection of such keV-scale candidates. The best studied models in this class are 
the sterile neutrino and gravitino, examples of fermionic superWIMPs [3l H], where it is 
well-known that there are no chances for direct detection. In contrast, the present paper 
studies bosonic superWIMP models, which need not entail such a pessimistic conclusion. In 
particular, we show that the potential sensitivity of direct detection experiments to bosonic 



keV dark matter can indeed be competitive with astrophysical bounds as well as with the 
lifetime and cosmic gamma-background constraints. 

The majority of underground direct-detection experiments are specifically tuned to the 
nuclear recoil energy regime and all 7- or /3-like events are typically considered as back- 
ground. Currently, the experiments with the best limits are now sensitive to picobarn-scale 
cross-sections for WIMP-nucleus elastic scattering, when the WIMP mass is on the order 
of the electroweak scale. To illustrate our main point, we present a simple estimate that 
compares the rate of absorption of 10 keV-mass axion-like dark matter (a-particles) , coupled 
to electrons with strength me/ fa = 171^(10^^ GeV)~^, with the recoil signal from a 1 TeV 
WIMP that has an elastic scattering cross section on the nucleus of cXez ~ 10~^^cm^: 

-....^DM-. . mwiMP ^ ^ ^ (^0'°f^\)" ^ io« X 10^ X 10-^^ ~ 0(1). (1) 



O^elVDMUwiMP "^a ^DM 10^^ Cm^ 

For this rather crude estimate we simply took (TabsVuM ~ fa"^, and maUa ~ ''t^wimp^wimp ~ 
Pdm- Although the superWIMP absorption cross section is orders of magnitude below the 
weak cross section, it is compensated by the tremendous gain in the local number density 
of dark matter particles, and by the fact that the inelastic cross section aabs scales inversely 
with the dark matter velocity. Although rather imprecise, the estimate ([1]) illustrates that 
the absorption of 10 keV dark matter can indeed produce signals that are well within modern 
detection capabilities. Moreover, the scale /„ ~ 10^° GeV is at or above the limits imposed 
by even the most stringent constraints on star cooling. We believe that this is an impor- 
tant point, and the sensitivity to new and viable dark matter scenarios can (and should) 
be explored by direct detection experiments. For comparison, fermionic candidates in the 
superWIMP class, a well-studied example of which is sterile neutrino dark matter, do not 
fall into the class covered by the estimate ([T]) for many different reasons; the primary one 
being that sterile neutrinos are not fully absorbed but rather converted to active neutrinos 
that carry away most of the rest energy. 

Thus far, the issue of direct detection of keV dark matter has only been discussed by the 
DAMA collaboration [5)^1 in connection with an annual modulation signal initially reported by 
DAMA/Nal [7j and recently confirmed by DAMA/Libra [8] . If indeed the annual modulation 
signal of DAMA were to be attributed to the absorption of keV dark matterp which would 
cause ionization but no significant recoil, it could provide a plausible explanation for why 
other front-running experiments such as CDMS and Xenon [111 [12] see no signal. This 
was the main point of Ref. [5], and in this paper we provide a re-analysis of this possibility, 
reaching instead a negative conclusion: the models that were presented in [5] as explaining the 
modulation are in fact ruled out either by lifetime arguments, by astrophysical constraints, 
or directly by the large unmodulated counting rates in underground detectors. 

To make the discussion sufficiently general, we focus on three generic possibilities for 
bosonic superWIMP dark matter: pseudoscalars, scalars and vectors, to be defined more 
precisely in Section 2. Some of these models are technically natural in the sense of hav- 
ing protection for the light dark matter mass, either by symmetry alone or by symmetries 



^As this paper was being prepared for publication, another experimental collaboration, CoGeNT [6], 
reported their results on the absorption of keV-scale superWIMPs. 

^See O [6] for recent analyses of WIMP models in relation to the DAMA results. 



combined with the imposed smallness of their couphng to the visible sector. For each of 
these models, in Section 3 we analyze the lifetime, emission, and absorption rates that de- 
termine the dark matter abundance, the level of diffuse and galactic gamma backgrounds, 
the efficiency of star cooling, and the rates for direct detection. The appropriate parameter 
space for each model is considered in Section 4, where we impose the relevant astrophysical 
constraints and determine the viability of direct detection. We also comment on the low 
level of any modulated component of the signal, and conclude with some additional remarks 
in Section 5. 



2 Light dark matter candidates 

In this section we list the dark matter candidate scenarios to be considered. 

• Pseudoscalar DM 

We start with the pseudoscalars a, and write the interactions as a combination of several 
derivative-like operators of dimension five: 

/:int = '^F.^F^'^ - ^^n'i^ + ■■■ (2) 

Ja Ja 

where F^,y and ip are the electromagnetic field strength and the Dirac field of the electrons, 
and the ellipsis denotes possible interactions with other fermions and gauge bosons, and 
for simplicity we shall assume a similar strength for the a-SM couplings in those sectors. 
Notice that the other possible pseudoscalar coupling to the electron, aipi'j^ip, can always be 
decomposed into the two operators in ([21) using the equations of motion once we account for 
the chiral anomaly. While the dimensionful coupling fa does regulate the overall strength of 
the SM-a interaction, the dimensionless coupling to photons C^ is crucial for determining the 
lifetime and 7-background created by a decay. Restricting ourselves to the electron-photon 
sector, we expect three generic possibilities for the size of the coupling C^: 
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Case A corresponds to a pseudoscalar coupled to photons at some UV normalization scale, 
with couplings to electrons generated radiatively. Having normalized the electron coupling to 
1/ fa in ([2]), the coupling C^ S> 1 in this case. Case C is the inverse of Case A. The derivative 
coupling to the electron axial-vector current may only lead to the F^^F^'^d'^a operator at 
loop level, hence the {m^a/rneY suppression. Finally, Case B is intermediate, when a is 
initially coupled to the fermion via the m.ijji'jQilj pseudoscalar operator. Clearly any of these 
three choices can be realized without fine tuning given an appropriate UV completion. An 
additional advantage of this model is the automatic protection of the pseudoscalar mass 



against radiative corrections, exactly as in the conventional axion case. As we are going to 
see later, only option C allows for the possibility of keV-scale dark matter without imposing 
overly strong constraints on the size of fa- 

• Scalar DM 

A similar looking Lagrangian can be written for the scalar case: 

/:int = ^F^^F'^-' - ^niei^^ + ■■■ (4) 

Is Is 

Here there is clearly no protection for the mass against radiative corrections. However, one 
can still exploit the smallness of the coupling f~^, to render a keV-scale mass technically 
natural. A one-loop correction will typically induce a mass term that scales as 

A(m.)= ~ ^. (5) 

where ruf is the mass of the heaviest fermion, and Auv is the ultraviolet cutoff. Taking both 
to the weak scale (implying supersymmetry) , and requiring rris ~ 10 keV is equivalent to 
having fg > 10^ GeV, in other words right at the boundary of the interesting regime for the 
couplings. As for the couplings to photons, both Cases A and B cases are plausible, while 
case C is tricky and requires some fine-tuned UV physics to cancel the main contribution 
from the rUe threshold. Finally, we note that the simplest renormalizable and SM-gauge- 
invariant realization of (jll) is to have the scalar singlet s coupled to the Higgs doublet via the 
relevant operator sH^H. We will not consider the scalar example in detail in what follows, 
but in many ways its phenomenology is similar to that of the pseudoscalar case. 

• Vector DM 

Finally, we introduce a model of keV-scale vector dark matter. We choose the intitial 
Lagrangian in the form identical to that studied in ^0\, where an extra U(l)' gauge field is 
coupled to the SM via kinetic mixing with the hypercharge field strength, 

1 K 

^ = ~~^tiv " -^ii.vF^t.v + Ch' + £dim>4, (6) 

with Chi encoding the physics responsible for breaking the U(l)', and £dim>4 includes possible 
non-renormalizable higher- dimension interaction terms such as if^ifV^^, y?lyF\^ etc. After 
the breaking of this secluded U(l)', the model takes the simplest possible form, 

^ = -^-^l. + \rn\yl + «:K9^F^., + ■ ■ ■ , (7) 

where we retained only relevant and marginal operators, and suppressed the U(l)' Higgs 
sector. This is one of the simplest UV-complete extensions of the SM, and it has been 
addressed in connection with electroweak-scale physics on a number of occasions [131 flO] . 
Most recently a study of cosmology in this model was performed for small sub-eV values of 
my [H| . For vanishingly small values of m^ , the extra sector decouples as mixing with the 





Figure 1: Dominant decays to photons, (a) 2-photon decay of the pseudoscalar a, and (b) the 3-photon 
decay of the vector V. 

photon can be reabsorbed into the mass term. This leads to an additional suppression of 
'j-V conversion at temperatures much in excess of mv [H], but for keV-scale dark matter 
this issue is often less important, and ([7]) can simply be traded for 

^ = -^V^u + \ri\yl + e«:KV'7MV^ + ■ ■ ■ , (8) 

where V^ couples to the electromagnetic current. For convenience, we introduce an analogue 
of the electromagnetic coupling and of the fine structure constant, 

e' = eK, a' = ^, (9) 

that necessarily appear in all rates for the emission or absorption of U(l)' vectors by SM 
particles. 

Obviously, the mass of the f/'(l) gauge boson is protected by gauge symmetry, and 
the question of naturalness is relegated to the corresponding U(l)' Higgs sector. However, 
given the smallness of the couplings that we are going to consider, e.g. k ~ O(10~^'^), the 
naturalness problem is no more severe than for the SM Higgs. 

The coupling of V to neutrinos is also possible via mixing with the Z-boson. This mixing, 
however, is further suppressed by a factor of {nia/MzY ~ 10"^"^, and we will disregard it in 
the analysis of the model. Other realizations of U(l)' [151 [IB], such as a gauged version of 
B — L, would allow couplings to neutrinos of the same size as the couplings to electrons. 



3 Decay, emission and absorption of superWIMPs 

3.1 Pseudoscalar (and scalar) DM 

Decays: Given an a-boson (or s-boson) mass below the electron threshold, the decays will 
be almost exclusively to photons, and mediated by the 2-photon interaction in ([2]), with the 
appropriate low energy value of /„. The decay is shown in Fig. 1 and the width is 

r-2. = ^m^. (10) 

Requiring that the dark matter lifetime be at least the age of the universe implies, 

r„r„^., < 1 ^ C^ < 2 X 10- X [J^y X (12i!X)'. (U) 



which we see is already a significant constraint for Case A in particular, and also for Case 
B if we stick with fiducial values for fa and ma- In constrast, Case C is less constrained. 
However, we also need to consider the 7-background created by a-decays and, as we will see 
in the next section, this provides a more stringent constraint (as is also the case for majoron 
models of dark matter, discussed recently in [T7]). 

Emission: An important process to consider is the a-emission from thermal states. This is 
relevant for determining the relic cosmological abundance, but is also the source of important 
constraints arising from new energy-loss mechanisms in stars. To obtain an estimate of the 
impact on star-cooling, we will focus on the Compton-like process e + 7 — *> e + a. Working 
in the limit ma,uj -C rrie, we obtain the following cross section. 



'^67— >ea '3^ 9 „r) ^ 






(12) 



where to is the photon frequency, and Va = (1 — ml/cu'^Y^^ is the velocity of the outgoing 
massive axion. In the limit ma — > 0, fa — *> 1, it reduces to a well-known result in conventional 
axion physics (See, e.g. [18] and references therein): 

3 m2/2 

The cross sections ( !T2l) and (TT3|) for pseudoscalar production can be translated to an 
energy-loss fiux (energy/volume/time) relevant for solar and red giant physics. In the limit 
of small axion mass we can estimate this rate as 

*_.„..(.._)^(^.')(|.)(^|^), .4) 

where pp is the electron Fermi momentum. For ma > T this formula needs to be sup- 
plemented by a factor of exp(— ma/T) to account for the Boltzmann-suppressed fraction of 
photons with energies above the axio-production threshold. We will compare this fiux to 
various constraints on stellar energetics in the next section. 

One of the most stringent astrophysical constraints on exotic particles often comes from 
supernova (SN) physics. Owing to the high temperature scale ~ 0(10 MeV) during the 
explosion, the coupling to electrons is not the dominant mechanism for axion production, as 
the rate effectively receives an additional suppression by a factor of m'l/T'^. It is well known 
that the coupling of a to nucleons provides far better sensitivity p^, and for the purpose of 
making an estimate we shall explore faqq ~ faee ~ fa- Since we consider pseudoscalars in the 
keV mass range, Tsn 3> m^a, and the emission of axionic dark matter particles in supernovae 
then differs little from the standard case of "invisible" axions [19j. 

The same argument applies to the thermal emission of axions in the early Universe. 
This occurs due to the interaction of SM fermions in the primordial plasma, and scales as 
Fq^ ~ Tm.'^f~}^ as long as the fermionic species ip is present, T > m^, and ip couples to the 
axion with a dimensionless coefficient ~ m,^/ fa^^- Taking into account that the Hubble rate 
scales as iJ ~ N^''^T'^Mp^, where N is the effective number of degrees of freedom, we can 



estimate the resulting number density of pseudoscalar particles weighted by the entropy, 

--1^ jnif-l^B,,^^^j^ (15) 

Vi "^ ip ■'"■vv 

where in the last relation we took into account that n^/s is maximized near the annihila- 
tion/decay threshold, T ~ m^, and the numerical constants B^ are introduced to account 
for specific details of each ip threshold. It was also assumed that pseudoscalar production 
happens with a rate slower than the Hubble expansion, so that the a-bosons are never fully 
thermal. Assuming no hierarchies among the couplings to the SM fermions, one observes that 
the largest contributions come from the heaviest fermions, which presumably will be the top 
and bottom quark. In principle, all B^ parameters can be calculated exactly, but we will not 
pursue this here as it will take us too far afield and the resulting formulae will still contain 
the model-dependent factors faqq- Instead, we will restrict ourselves to a simple dimensional 
estimate of the dark matter axion abundance produced via the 6-quark (production from the 
top quark depends on other details such as the electroweak phase transition): 



X — X — ~ 5i?5 X -— — X . (16) 






f^baryon 1 GeV S Ub keV V fabb J 

Since Bf, is naturally of 0(0.1 — 1), one can see that thermally generated peudoscalar dark 
matter has an abundance in the right ballpark given a keV-scale mass and the most inter- 
esting range for the couphng, fabb ~ 10^° GeV. A precisely analogous argument holds for 
the scalar super- WIMP abundance, where a coupling scale in excess of 10^° GeV is further 
supported by the technical naturalness argument, Eq. (^. 

Absorption: For the energy range considered here, the axioelectric variant of the photoelectric 
effect is the most important process. A calculation of the axioelectric effect on an atom was 
performed earlier [20] in connection with the possibility of detecting axions with keV energy 
emitted by the solar interior. We would like to rederive the same effect in our setting, 
specializing to absorption of the nonrelativistic keV mass pseudoscalar. Working to leading 
order in Va/c, one can reduce the interaction in ([2]) to the following term in the nonrelativistic 
electron Hamiltonian: 

£,« = -^*7''7V ^ i/,« = ^^, (17) 

fa fa rUe 

where p is the momentum operator for the electron. Using a{t) ~ exp(— -irriat) and the 
nonrelativistic Hamiltonian, the matrix element of ifjnt reduces to the following expression, 

9 
HI 

Mf, = -^{f\{rcT)\t), (18) 

Ja 

where i and / are initial and final state electron wave functions with Ef — Ei = rria. This is 
analogous to the amplitude for the El absorption of a photon, M^ = uj{f\e{ve)\i), with the 
photon polarization e exchanged for the spin operator. Thus the axioelectric effect is very 
similar to the photoelectric effect with photon energy a; = m^. The only difference is that 



the wave function of the absorbed photon contains a space-dependent factor exp(ikr) with 
k = uj, which in the absorption of a massive vector particle is replaced by exp(imavr). Here v 
is the velocity of the incoming DM particle. Since v ~ 10~^c, the latter oscillating factor can 
safely be taken to 1. The effect of the photon spatial momentum is however parametrically 
small in comparison with the momentum transfer to the final electron, 0{y/m^Tn^), as long 
as the energy in these processes is much smaller than the electron mass mf.. Neglecting the 
minor effect of the photon spatial momentum, squaring flTSl) . averaging over the initial spin, 
and summing over the final, we arrive at the following approximate relation between the 
cross sections for axion absorption and the photoelectric effect: 



(^absV 3 ml 



(19) 



o-photoi^ = ma)c 4:71 a f^ 

Notice that ratio (TT9l) remains finite even in the limit of f -^ 0. 

For the purposes of an experimental search for pseudoscalar dark matter, it is useful to 
express flT9|) directly as a counting rate in a detector consisting of a single atomic species of 
atomic mass A, 

1.2 X 10^9 fma\ faphoto\ , _i , _i .„„. 

where we used the local dark matter density, p^M = 0.3 GeV cm^^, and introduced the 
dimensionless coupling Qaee = 2r7ie//a to allow for direct comparison with the existing results 
of DAMA [5l E] and the most recent experimental paper by the CoGeNT colaboration [6] . 

Note that our results differ from those quoted in [5] and [6]c| We observe that the 
expression used in |5j for the axio- ionization cross section results from the omission of the 
leading term in the axion-electron Hamiltonian (TT71) . Indeed, Eq. (45) of Ref. [5| has an 
interaction term that vanishes as vdm —>■ giving a subleading contribution to the counting 
rate suppressed by f|)^j ~ 10^^. We shall return to this issue in section 4.2. 

For completeness, we will also quote the axioelectric cross section for massless axions 
with energy of O(keV), as this result is in widespead use for the solar physics of axions. 
Taking Ua ^ rrie, and rria — i> 0, we arrive at the following analogue of the matrix element 

(HHD, 

M/. = ^(/|(r^)-(n^)(nr)|.), (21) 

Ja 

where n is the direction of the incoming axion. Following the same route as before, we arrive 
at a relation between the photo- and axio-electric cross sections: 



O'abs 



O'photoi^ — ^a) 



'^'^ (22) 
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We observe that this result is twice as large as the formula regularly quoted in the literature 
(see e.g. [20]). The source of the discrepancy with previous calculations can be traced to the 



^Following the release of version 1 of the present paper, the rate formula used in [6] - which originally 
made use of the relativistic expression for aabs of [20] - was corrected in 0807.0879v4 to use the appropriate 
non-relativistic expression (fT9|) . with a result for R consistent with ([20]). 



absence of the first term in Eq. (pTj) . and this corrected formula may prove useful for solar 
axion searches. 

Finally, on a more pedagogical note, we would like to demonstrate explicitly that the 
alternative choice for the pseudoscalar coupling to electrons, {2me/ fa)ciipi'y^ip in £int, leads 
to the same expression for Hi,^t and the matrix elements, as the axial vector coupling (IT7|) . 
Since both forms are related up to the total derivative, this must necessarily be the case, 
but it is useful to see how this works in detail. We will only consider the case of a massive 
nonrelativistic pseudoscalar as is relevant for this paper and, using the v/c expansion, we can 
then write the lower component x(r) of the Dirac spinor ip in terms of the upper component 
0(r), the total energy E, the potential energy f/(r) and the helicity operator (per), 

2mex(r) = (^1 - ^^^^ (P'T)0(r). (23) 

With this expression, the matrix element takes the form 

M, ^ ^««..,), ^ l(/|,p.) (l - ^^) - (: - ^^) ,p.,|,. ,24) 

Given the Hamiltonian Hq = p^/(2me) + U{r), with HQ\i{f)) = Ei(^f)\i{f)) , straightforward 
quantum mechanical manipulations reduce the matrix element to 

= -£^(/|(p^)K) = ^(/l(rtT)K), (25) 

where we have used Ej — Ei = rria. The result is identical to (fT8|) as expected. 

3.2 Vector DM 

Decays: The vector model of keV dark matter has an important distinction when compared 
to scalar/pseudoscalar models, as the direct decay to two photons is strictly forbidden, 
regardless of how the new vector particle is introduced in the model. 

The decay to three 7 quanta is allowed at the loop level. In the limit my -C rrie, the 
electron coupling will then generate a dimension-eight interaction with photons of Euler- 
Heisenberg form, 

3 / 
^7 = ^^^^2^4 {UF^uF,,F,,V^^ - 5{F,,F^,){F,,Vx^)) , (26) 

where F^^, is the electromagnetic field strength tensor and V^^, = d^Vy — d^V^ is its analog 
for the field V^ of the massive vector. These operators then mediate the dominant 87 decay 
below the electron threshold, shown in Fig. 1, and we will calculate the spectrum and decay 



rate explicitly. We find that the full Dalitz plot distribution of photon energies in the rest 
frame of V is described by 

-556 [ujI{u02 + ^3) + OjI{uOi + CJ3) + U0l{u0i + UJ2) - UJiUJ20J^{l^l +UJ2 + ^3)] } J (27) 

where cJi, uj2 and uj^ = my — uJi — UJ2 are the final-state photon energies. Integrating over 
one of the photons yields the inclusive one-photon spectrum in the decay, 

dr a^a' ml , / 2919 



dx 27375%3 ml 



( 2919 A 

( 1715-3105X+— -xM , (2^ 



where x = 2uj/m,y, so that the physical region for the dimensionless parameter x ranges 
from to 1. Integrating the one-photon spectrum we finally obtain the total decay rate: 

17 aV < ^ (4 70 X 10-«) a'a' ^ f29) 

An immediate application of this result is to provide a constraint on the masses and 
couplings such that the dark matter hfetime is at least equal to the age of the universe. We 
find 

TvTv^3^,<l =^ my(a')'/'<lkeV. (30) 

Even if the lifetime constraint is satisfied, the 7 background created by V decays can be 
detected as a diffuse cosmological 7 background, or as an extra contribution to the 7 back- 
ground in our galaxy. 

Emission: In considering emission, we will again focus on Compton-like scattering with 

electrons, e + 7 -^ e + V, which will be relevant in considering energy loss in stars. In 

the limit where m,v is small compared to the frequency of the absorbed photon, and in 

the nonrelativistic approximation with respect to the electron, this process has a standard 

Thomson-like cross-section: 

Siraa' 

3m^ 



OIIU.U. 
CTej-^eV — — ^- [oi-) 



This results in an energy loss flux (energy/volume/time) from a thermal plasma which we 
can estimate as 

to again be compared with various constraints on stellar energetics in the next section. 
As above, this rate has to be modified to account for Boltzmann suppression if m,v ^ T. 
However, the formulae (13T|) and (132|) neglect another possibly important effect: the effective 
suppression of V^ — 7 mixing due to the dynamical mass of the photon inside plasma. This 
screening means that inside a highly-conducting medium the coupling constant a' has to be 
modified as follows: 

m,Y \ . a m,y 

— 5 ^ -^ Opff ~ — 7— lor iTT'V < mn- (33) 

rriy — mjj J m^ 
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For temperatures which are low compared to the electron mass, the effective dynamical mass 
(or plasma frequency) can be taken as m'jj = A'Kam~^{nf. + ne+). In the stellar environments 
where ne+ = 0, and rig is not exceedingly large, this would not imply a strong suppression of 
the emission rates unless my is under a few keV. However, for cosmological or SN applica- 
tions, where T > irte-, the high-temperature effective mass is relevant, m^ = AttoT"^ /?>. The 
consequent suppression of all \^-production rates can then be quite significant. 

To estimate V^-production during supernovae, we note that it is the coupling of V to 
nucleons, and in particular the coupling to the neutron magnetic moment that leads to the 
most important production channel. This logic is motivated by the large number density of 
nuclear matter in the core of the SN, and by the analogy with the axion case, where coupling 
to nucleons provide more stringent constraints. Furthemore, since the nuclear matter in the 
core is mostly neutrons, it is the neutron electromagnetic formfactors, and magnetic moment 
in particular, that mediate \^-pro duct ion. Adopting the method of Refs. [2T[[22], we estimate 
the emission of vectors by factorizing the nucleon elastic scattering cross section a^N and 
the probability of ^-emission due to the neutron spin flip, 

^NN^NNV ^-^^(^NNnlf ( ) Ta'^^{T), (34) 

47r^ \mN J 

where /xat = 1.9(47ra)-'^/^/(2mAr) is the magnetic moment of the neutron, and njq is the 
neutron number density. Assuming typical temperatures on the order of 10 MeV, we notice 
that for my < 1 MeV, the m^i in (!33l) dominates, and it is possible to reduce (IMI) to 

^NN^NNV- 10 X - X ^5/2^,;, > (35) 

which has a strong dependence on my but a rather mild temperature-dependence. One 
observes that the primary advantage of the SN limits, namely an enhancement of the emission 
rates at large temperatures, is completely lost in ( l35l) due to the strong suppression of a'^g. 
Consequently, we do not find any competitive constraints on the parameter space from SN 
physicsjj The inclusion of other V^-production channels {e.g. from electrons, positrons, and 
other charge carriers) would not change this conclusion. 

The cosmological abundance of keV-scale vectors has a thermal and also possibly a non- 
thermal component. Thermal emission of vectors due to the operator nF^uVfj^,^ is strongly 
inhibited at high temperatures, and thus the most important threshold for emission is the 
e~^e~ threshold at T ~ O.SMeV. Production of V occurs via the processes e + •y ^ e + V 
and e + e"*" — > y + 7, and even below the the electron threshold via j + 'j —>■ 'j + V mediated 
by a loop diagram, although the latter process is probably going to be subdominant. In this 
paper we will concentrate on the mass range my < 100 keV. This allows for some additional 
simplifications, as the production rate is peaked below T = m^ due to (133|) . For my = 10 
keV, the peak in production occurs at ~ 150 keV, which is predominantly due to Thomson- 
like scattering off the remaining electrons and positrons. This allows for a nonrelativistic 
treatment of the electrons, and justifies the use of the Thomson-like formula. With these 



* Additional suppression of ([55]) arises from the large residual chemical potential for electrons in SN 
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simplifications, we can calculate the freeze-out abundance of V^-particles per photon as a 
function of niv and a', 

Yv{mv,a) = — ^ —2 ^ x — — / , . , Jo,.^^ , (36) 



n. 
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HT a "iml Jq {m\ + mlYHT 



where we took into account that the number densities of electrons and positrons are approx- 
imately equal, Ue — ne+ — 2"^/^(meT/7r)^/^exp(— me/T). Once again, we assume that the 
^-sector is not in chemical equilibrium with the SM, and the production rate of ^-particles 
is slower than the Hubble expansion rate H. The upper limit of integration is chosen as 
~ rrte, the scale where the nonrelativistic approximation breaks down, but can formally be 
extended to infinity because the integrand has a maximum below T = rrie and falls quickly at 
high temperatures. Using (1361) . we can immediately calculate the corresponding contribution 
of V-particles to the total energy density today: 

nvimv, a') = fibaryon — ^^^^ Yy = 73 Yy X ^. (37) 

rriN nbaryon kcV 

Choosing mv = 10 keV for example and tuning Qy to the measured value of ^dm, 
i.e. (]y(10 keV, a') = 0.2, we find 

- ~ 10-20. (38) 

a 

This tiny coupling certainly justifies referring to V dark matter as a super- WIMP. 

Eqs. ( l36l) and ( 1371) provide a way of estimating the thermal component of V dark mat- 
ter created at the electron threshold. There are, of course, other possiblities for creating 
additional contributions to Qy Foi' example, infiation may end with some infiaton decays 
to particles in the U(l)' sector; higher dimensional operators may provide an efficient way 
of transfering energy from the SM sector into the U(l)' [21]; the Higgs'-strahlung processes, 
e"*" + e~ — ^ V"* — *> V^ + h', might also be important, etc. All of these mechanisms are rather 
model-dependent, and can only add to fly on top of the estimates (15B]) and (P7|) . 

It is worth dwelling on the Higgs'-strahlung process for V^-production, as it can be an 
important mechanism for the following reason. The virtuality of the 'y — V line is g^ > 4mg, 
and thus one does not have any thermal suppression for this process. Moreover, the cross 
section for Higgs'-strahlung remains finite in the my -^ limit, and in our model scales as 
~ a'(e)2/i?2, where E is the center of mass energy, and e is the gauge coupling in the U(l)' 
sector. For a fixed value of e, the smaller the mass of my we consider, the more important 
the Higgs'-strahlung production of V will be. In this paper, we choose to ignore it noting 
that for our choice of mass range, my > keV, the thermal suppression fl33|) at the electron 
threshold is at most a factor of 0(10), and we can suppress the Higgs'-strahlung relative to 
the Thomson production of V by choosing e < eo 



^ We note that in the recent paper [M], the cosmological abundance of V bosons for mv < 1 eV is 
computed, but without including the Higgs'-strahlung process which is likely to be very significant for this 
mass range. Indeed, Ref. J.4J often considers rather large mixing angles k, which would completely populate 
V and h' right at the electron threshold, creating a minimum of four new degrees of freedom, unless ne is 
chosen to be less than ^^ lO^^e. Having four new degrees of freedom at T = 1 MeV is now firmly excluded 
by Big Bang Nucleosynthesis constraints. 
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Absorption: The absorption of V DM by atoms is very similar to the photoelectric effect 
with photon energy u = niy. The only difference is again the factor exp(ikr) incorporating 
the photon spatial momentum. As discussed previously in connection with the axioelectric 
effect this factor can be safely approximated by one. Thus, for our estimates it suffices to 
take 

r-^-" , ^ -. (39) 

(^photoi^ = mv)c a 

Converting this cross section into a counting rate gives 



(40) 



^ ^ 4X10- <y /keV\ /^-j ^^-^A,y-K 
A a \ mv / V bn / 

4 Direct detection of keV superWIMPs 

In this section, we use the various interaction rates determined above to assess the viability 
of direct detection of bosonic keV-scale superWIMP dark matter. We will focus on the 
photoelectric-type ionization cross-section as the primary source of a direct detection signal, 
and also consider how various indirect constraints will cut into the available parameter space. 
In this regard, we consider energy loss from stars and supernovae, and also the 7-background 
produced by decay. We will now consider the pseudoscalar and vector models in turn. 

4.1 Pseudoscalar DM 

As has already been noted above, the monochromatic decay a -^ 77 leads to a rather 
stringent constraint on Cy/ fa not just through the need to have a sufficiently long lifetime, 
but more significantly through the induced source of galactic (and also diffuse cosmic) X- 
rays. The constraints on such monochromatic sources in the galaxy are quite stringent, and 
for this reason we will adopt the most conservative model. Case C, in relating C-y to the scales 
fa and ma- This allows us to restrict our attention to a 2-dimensional parameter space, fa 
vs rUa, which we show in Fig. 2. The various contours on this plot are described below: 

1. The monochromatic 7-decay leads to a very strong constraint from observations of the 
galactic background, and in particular searches for various line sources in this energy 
range [25]. For example, the decay width to monochromatic photons in the 10 keV 
range is constrained to be less than approximately F < 10~^^s~^, which we observe is 
around ten orders of magnitude more stringent than the constraint on the lifetime. As 
noted above, the strength of this constraint means that we will only consider Case C 
in which the photon coupling is suppressed. The limits obtained in [25] can then be 
used to form an exclusion contour for fa in terms of rUa as shown in Fig. 2. The growth 
of the decay rate with mass is sufficient to render the constraint more stringent at the 
upper end of the mass range considered. The couphng constants for Cases A and B as 
well as for scalar dark matter are even more strongly constrained. 
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Figure 2: Fixing C^ in terms of fa and rua according to Case C, we plot the direct detection sensitivity 
to pseudoscalar DM arising from the axioelectric cross-section on Ge, assuming a fiducial sensitivity of 
the detector equivalent to a Ipb cross-section for a 100 GeV WIMP. We also show the constraints arising 
from the He-burning lifetime in HB stars, from SN cooling via a coupling to the neutron magnetic moment 
with faNN = fa, and most significantly the monochromatic 7-background from decays in the Galaxy. The 
grey shaded region is excluded by the latter indirect constraints. The thick black line corresponds to the 
parameters required to reproduce the required dark matter abundance from thermal production with fabb = 

Ja- 



2. An important constraint on any new light states that couple to photons or electrons is 
that the new energy loss mechanisms should not severely disrupt the life-cycles of stars. 
Earlier, we used a Compton-like emission process to estimate the energy flux into a 
particles. We will consider two constraints, the first of which constrains energy loss 
from He-burning Horizontal branch (HB) stars in globular clusters. Following Raffelt 
|19] (see also Raffelt and Weiss ^^ for constraints from He ignition), this constrains 
the energy flux to $ < 10~^^ MeV^ at a density of p ~ lO^g/cm^ and a temperature 
of T ~ 10 keV. We plot the corresponding exclusion contour in Fig. 2, which degrades 
for rria ^ T due to the Boltzmann suppression of photons with uj > irta. 

An energy loss constraint which is less sensitive to the mass, at least in the keV-range, 
arises from cooling of the supernova core, e.g. in SN1987A. In this case, the constraint 
on the flux is $ < 10~^^ MeV^ at a density of p ~ lO^^'g/cm^ and a temperature of 
T ~ 30 MeV [inj. At this core temperature, the a particles are effectively massless 
and from the electron coupling we find a rather weak horizontal exclusion contour of 
fa > 8 X 10^ GeV. However, allowing for derivative couplings to quarks, in analogy 
to the electron coupling considered above, a stronger constraint on fa ~ faNN ensues 
from the induced coupling to the neutron magnetic moment in the degenerate SN core 
leading to the contour shown in Fig. 2. 
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3. We also include a line that corresponds to Qa = ^dm, using Eq. (fT6l) . assuming 
Bh ~ 0.1 and choosing fabb = fa- This "natural abundance" hne competes with the 
SN constraints, but is already deeply inside the area excluded by the galactic X-ray 
constraints. 

4. The axioelectric sensitivity contour assumes that an experiment like CDMS has a 
sensitivity to the ionization signal from absorption which is equivalent to its sensitivity 
to the recoil of a 100 GeV WIMP with a cross-section per nucleus of Ipb. This is 
simply a benchmark point, and the contour can be rescaled according to the fact that 
(av) oc /~^. This sensitivity line is derived under the assumption that the a-bosons are 
the dominant component of the dark matter energy density in the solar neighborhood, 
which again can easily be rescaled to a more generic case. Note that we have made 
use of the well-measured photoelectric cross-section on Ge, which has a sharp break at 
around 11 keV, which is smoothed out somewhat in the contour. 

4.2 Vector DM 

For vector DM, gauge invariance is rather important in restricting the photon decay rate, and 
we observe in this model that the indirect constraints particularly from the 7-background 
are considerably weaker. This model has only two parameters, and we plot the space a' /a 
vs mv in Fig. 3. The various contours are the same as those described for pseudoscalar DM, 
and use the same constraints as above. We briefly describe the distinctions below: 

1. The galactic gamma background is less dominant in this case, and the constraint here 
is an estimate that comes from relaxing the bound on monochromatic lines by an order 
of magnitude to account for the broader, but still quite peaked, distribution for the 
37-decay. The constraint is again stronger for larger mass. 

2. The stellar bounds again arise from energy loss due to the Compton-like processes dis- 
cussed in the previous section. We observe that the stellar constraints are particularly 
strong for TUy of order the core temperature of HB starqj, but degrades exponentially 
for my > 10 keV due to Boltzmann suppression, while in this case SN physics does 
not provide strong constraints anywhere on this plot on account of the effective photon 
mass suppression of a^g. 

3. The V-electric sensitivity was obtained in the same way as for the pseudoscalar, and 
we see that in this case the prospects for direct detection look particularly strong, 
and indeed this approach may have the best sensitivity by an order of magnitude if 
the relevant experiments can discriminate the ionization signal. As before, this line 
is calculated assuming that V bosons are the dominant component of galactic dark 
matter. 



^We thank G. RafFelt for emphasizing the importance of constraints from the He-burning hfetime of HB 
stars in these models. 
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Figure 3: We plot the direct detection sensitivity to vector DM from the V-electric cross-section on Ge, 
assuming a fiducial sensitivity of the detector equivalent to a Ipb cross-section for a 100 GeV WIMP. We 
also show the constraints from the He-burning lifetime in HB stars, and the 7-background from 37-decays 
in the Galaxy. The grey shaded region is again excluded by the indirect constraints, while the thick black 
line corresponds to the parameters required to reproduce the required dark matter abundance from thermal 
production. 



4. Remarkably, part of the natural abundance line, calculated using (136!) and (!37|) . falls 
within the region allowed by the astrophysics constraints. It can nonetheless be probed 
rather effectively with direct dark matter searches. 



4.3 Discussion of annual modulation 

An annual modulation of the counting rate is a welcome feature in the direct detection of 
WIMPs. In the case of elastic WIMP-nucleus scattering, the cross section typically has an 
s-wave component, which is a constant independent of velocity. The counting rate, however, 
is proportional to the total dark matter flux and thus acquires a seasonal modulation due to 
the Earth's motion at the level of Avdm/vdm ~ 0.05. The DAMA/Nal and DAMA/Libra 
collaborations have utilized this idea to search for WIMP-nucleus recoils, using the annual 
modulation as a filter to separate signal from background. 



4.3.1 Suppressed modulation of the super- WIMP counting rate 

In principle, the absorption of bosonic superWIMPs could also generate an annual modula- 
tion of the signal [5|. This idea, however, immediately runs aground because inelastic cross 
sections are typically inversely proportional to the velocity of the incoming nonrelativistic 
particles [2S!- In combination with the flux, this renders the counting rate independent of 
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the velocity, and therefore un-modulated by the Earth's motion at the experimentally rele- 
vant percent level. We find that this is indeed the case for both classes of models considered 
in this paper, pseudoscalar and vector superWIMP dark matter, with rates given in fl2UI) 
and (HOl) . We can nonetheless estimate the degree of modulation in the absorption, noting 
that the dark matter velocity enters via exp(ikr), which in the Born approximation can 
be combined with exp(ipr) representing the wavefunction of the outgoing photo-electron. 
Therefore, modulation of the counting rate arises primarily due to the modulation of p. 
Since k ^p and the cross section is generally a smooth function of energy, we estimate that 

Modulated absorption ~ ~ — - ~ x \ma(v)nne- (41) 

p ^mema(v)C c V ' ' 

For the characteristic energy range of the DAMA signal, rria^v) ~ 5 keV, the modulation does 
not exceed 10~^. One could potentially worry that ma(v) Kiay turn out to be exactly equal 
to some ionization threshold where the assumption of a smooth cross section as a function 
of energy breaks down, and the momentum of the photoelectron is suppressed relative to 
its natural value. In this case the modulated part of the cross section can be enhanced, 
but not by the four orders of magnitude required to explain the DAMA signal. This is 
because the detector comprises many-electron atoms and unmodulated absorption by other 
shells would become important. From Eq. pTj) . one immediately concludes that the reported 
DAMA modulation signal cannot be explained by the absorption of axion-like or vector-like 
particles. Specifically, if the couplings are tuned in such a way that the modulated part of 
the axio-electric effect matches the DAMA signal, the unmodulated part will be in excess of 
the total number of events that DAMA observes by four orders of magnitude^ 

4.3.2 Solar neutrino backgrounds 

Leaving aside the issue of (un)modulated absorption of superWIMPs, we would like to dis- 
cuss other possible sources of annual modulation in ionization. For example, is it possible 
that emission from the Sun creates a modulated ionization signal? Although the photons 
themselves do not reach underground facilities, solar neutrinos can do so and create some 
amount of ionization. The size of the modulation at the level of ~ 0.03 — 0.05 [8] is roughly 
consistent with that of the solar neutrino fiux due to the eccentricity of the Earth's orbit. 
However, this modulation is expected to be in anti-phase with the DAMA result because in 
the (northern hemisphere) summer the Earth is farther away from the Sun than in the win- 
ter, and consequently the fiux of neutrinos is slightly lower during the summer months. Since 
solar neutrinos have energies much in excess of atomic ionization thresholds, one also does 
not expect a concentration of the ionization signal around a few keV, modulo unknown solid 
state effects. We can estimate the importance of this (modulated) neutrino background more 
generally. The contribution of neutral currents to elastic scattering on nuclei that mimic re- 
coil events has been estimated previously in [27j. It is only the most energetic fraction of 



^In the analysis of the axioelectric process in [5], we should note that, besides the omission of the leading 
term (I17|) in the axion-electron Hamiltonian, the subleading term '^ crk is also incorrectly averaged over 
the wavefunctions of the initial and final electrons. Since the axion momentum k is an external vector, the 
matrix element of the electron spin operator between two states with different energies can give a non-zero 
result only on taking the spin-orbit interaction into account, which is not consistent with [5]. 
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the solar neutrino flux that is capable of creating nuclear recoils with more than a 1 keV 
energy release. The rate of such events for DAMA should not exceed 10~^ kg~^day~^ [2Zj- 
Ionization may also be created by the main part of the pp neutrino flux. Our estimates 
show that the total counting rate due to ionization by pp neutrinos is also at the level of 
10~^ kg~^day~^, providing a modulated counting rate that is about two-to-three orders of 
magnitude smaller than the reported DAMA signal. Thus solar neutrinos do not induce a 
modulated signal capable of matching the amplitude observed by DAMA. However, for fu- 
ture high-sensitivity searches for bosonic superWIMPs, of the kind advocated in this paper, 
the ionization created by solar neutrinos may constitute an important source of background. 

To conclude this section, in a more speculative vein one cannot help noticing that the 
energy range for the modulated DAMA signal, 2 — 6 keV, claimed in [8] is comparable to the 
temperature of the solar core, and indeed matches the energy range where the emission of 
exotic massless particles would naturally be peaked [19]. This could be "standard" axions, 
or massless U(l)' bosons coupled to the SM via marginal operators [28], or another similar 
type of exotic. Again, the ionization signal created by these exotic particles can have a 
3% annually modulated component, but it will be "vr-shifted" relative to the DAMA signal. 
Only if the absorption within the Earth is somehow an important effect could the integrated 
day-night effect potentially induce modulation with a maximum in June and a minimum in 
December. We believe that this latter explanation can be directly checked using the DAMA 
and DAMA/Libra datasets. 

5 Concluding remarks 

With vast resources now justifiably being devoted toward the direct detection of dark matter, 
it appears all but clear that the scientific scope of these searches should be diversified. While 
a characteristic elastic WIMP-nucleus recoil may remain the main benchmark scenario for 
these experiments, the detection possibilities for other generic classes of dark matter should 
certainly be exploited. In this paper, we have shown that bosonic superWIMPs represent a 
legitimate and feasible target. 

The models analyzed in this paper can naturally produce the required relic abundance 
of dark matter, once the couplings a' or {Trtf^/ faY are in the superweak ~ 0(10"^'') range. 
The feeble nature of this coupling to the SM is partially overcome by a factor of (c/vdm) x 
{fnwiMp/fna{v)) which enhances superWIMP absorption relative to WIMP scattering and 
renders direct detection feasible. We have analyzed models that do not require any special 
fine-tunings and enjoy protection for the small mass scales due to the symmetries of their 
interactions: gauge symmetry for the model of secluded U(l)', and shift symmetry for the 
pseudoscalar model. 

The result of our analysis has revealed that for the pseudoscalar model, direct detection 
sensitivity may compete with the red giant, solar and SN constraints. However, we found 
that the model is most strongly constrained by limits on monochromatic X-ray lines in our 
galaxy, that rule out much of the interesting regions of the parameter space. In vector models 
of dark matter, the production of X-rays is strongly inhibited by gauge invariance, and for 
the model studied here direct searches for ionization are apparently capable of probing the 
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most interesting range of the coupling-mass parameter space, namely that consistent with 
the observed dark matter energy density. Interestingly, part of this range is not strongly 
constrained by indirect astrophysical probes, away from the low mass region where stellar 
energy loss constraints become important and the higher mass range where the 7-background 
is too large. Finally, contrary to some existing claims in the literature, we have found that 
the absorption of superWIMPs does not lead to an annual modulation of the ionization 
signal at a level that would be of experimental interest. 
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